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Thermal leptogenesis explains the observed matter-antimatter asymmetry of the universe in terms 
of neutrino masses, consistent with neutrino oscillation experiments. We present a full quantum me- 
chanical calculation of the generated lepton asymmetry based on Kadanoff-Baym equations. Origin 
of the asymmetry is the departure of the statistical propagator of the heavy Majorana neutrino 
from the equilibrium propagator, together with CP violating couplings. The lepton asymmetry is 
calculated directly in terms of Green's functions without referring to 'number densities'. A detailed 
comparison with Boltzmann equations shows that conventional leptogenesis calculations have an 
uncertainty of at least one order of magnitude. Particularly important is the inclusion of thermal 
damping rates in the full quantum mechanical calculation. 



Most theories of baryogenesis involve quantum inter- 
ferences in a thermal bath in a crucial manner Here 
we consider thermal leptogenesis [5| which, in its simplest 
version, is dominated by the CP violating interactions of 
the lightest of the heavy Majorana neutrinos, the seesaw 
partners of the ordinary neutrinos. For neutrino masses 
inferred from neutrino oscillations, leptogenesis is dom- 
inated just by decays and inverse decays of the heavy 
neutrinos in the thermal plasma 

Almost all leptogenesis calculations are based on Boltz- 
mann equations. This treatment has a basic conceptual 
problem: the Boltzmann equations are classical equa- 
tions for the time evolution of phase space distribution 
functions; the involved collision terms, however, are usu- 
ally zero-temperature S-matrix elements which involve 
quantum interferences. Clearly, a full quantum mechan- 
ical treatment is necessary to understand the range of 
validity of the Boltzmann equations and to determine 
the size of corrections 0, Q . 

Various thermal corrections @ have been incorporated 
in the context of Boltzmann equations, and 'quantum 
Boltzmann equations' have been derived from Kadanoff- 
Baym equations [tI, In 01, a solution of the Kadanoff- 
Baym equations for leptogenesis has been found to lead- 
ing order in a derivative expansion in terms of distribu- 
tion functions satisfying the Boltzmann equations. 

In this Letter we discuss leptogenesis directly in terms 
of Green's functions which are solutions of the Kadanoff- 
Baym equations, thus avoiding all approximations nec- 
essary to arrive at Boltzmann equations. Our work is 
based on |^9.], where the approach to thermal equilibrium 
has been discussed in terms of Green's functions for a 
toy model, a scalar field coupled to a large thermal bath. 
Here we extend this method to leptogenesis, which yields 
the lepton asymmetry directly in terms of Green's func- 
tions. In the following we describe the main result of our 
work. Detailed derivations will be given in [loj . 

The interactions of N , the lightest of the heavy Ma- 
jorana neutrinos, with the Higgs doublet <f> and lepton 



doublets I Li is described by the lagrangian (cf. [J), 



C = lu^Xl^N + KiClu^ - -MN' CN 

1 1 _ ~ _T ~ 

+ 2'^i]>'Li'l^ CIlj(I) + -jr]*jlLt(l) Cl^jcl) ; (1) 



here C is the charge conjugation matrix, 
the coupling 



fc>i 



— io'2(t>, and 



(2) 



is obtained by integrating out the heavy Majorana neu- 
trinos Nk>i with Mkyi 3> Ml = M . We shall consider 
the case of small Yukawa couplings, Aa ^ 1, such that 
the decay width of N is much smaller than its mass. The 
Lagrangian ([T]) represents an effective low-energy theory, 
valid for momenta up to Mk>i. 

The Boltzmann equations for the time evolution of 
the distribution functions of heavy neutrinos, Higgs and 
lepton doublets are well known In this Letter our 

main goal is the comparison of Boltzmann and Kadanoff- 
Baym equations. We therefore focus on the CP- violating 
source term for the asymmetry and ignore the washout 
terms and the Hubble expansion, which can be added in 
a straightforward way [10| . 

For the distribution function of the heavy neutrinos 
one has UM. 



d_ 

dt 



lN{t,Uj) 



/k,q 

X [/^(t,a;)(l-/Kfc))(l + /0(9)) 
-fi{k)f4q){l-fN{t,uj))], (3) 



where oj = \J + p^, fc and q are the energies of N , 
I and (/) with equilibrium distribution functions /; and 
Z^, respectively; the averaged decay matrix element is 
\M{N{p) ^ l{k)(j){q))\'^ =2(AU)^^p-fc (cf. For the 



2 



momentum integrations we use the notation 



L-I- 



(4) 



The sum of decay and inverse decay widths, which de- 
termines the rate for the approach to equihbrium jl3l] , is 
given by 

rp{u) = {x^x)^^- f {2i,)H\k + q-p) 

^ "'k,q 

P ■ k fi4,{k,q) , (5) 

where we have introduced the function 

Mk,q) = fl{k)Mq) + (1 - Mk)){l + U{q)) 

= 1- fl{k) + fM ■ (6) 

For the solution of the Boltzmann equation ([3]) with vac- 
uum initial condition, fiq{0,uj) — 0, one easily obtains 

(7) 

where /^'(w) = l/(e'^'^ + 1), and P = l/T is the inverse 
temperature. 

The Boltzmann equation for the lepton distribution 
function is given by 



{2^Y5\k + q~p) 



d 1 

a^^'(*'^) = -2fc, 

X [\Milcj, ^ N)\'fi{k)fM(^ - .fN{t,oj)) (8) 
- \MiN ~> lcj))\^fN{t,io){l - fi{k)){l + Uiq))] , 

where now ©(A^) corrections to the matrix elements have 
to be kept. Using Eq. ([7]) one obtains for the lepton asym- 
metry fLi = fii - fii, with initial condition /Li(0, k) = 0, 

fu{t,k) = -eul f {27T)H\k + q - p) p ■ k 



x/z0(fc,g)/^''H-(l-e-i^*) , (9) 



where 



3Im{A*i(7^A*),-i}M 



(10) 



Summing over all flavours, the generated asymmetry 
is proportional to the familiar CP-asymmetry: e — 
j:,eu/ (X^X),, = 3Im(At7yA)M/(167r(AtA)^^) fl. 

For later comparison with solutions of the Kadanoff- 
Baym equations, it is convenient to rewrite ^ in the 
form 



fu{t,k) 



167r 



k ■ k' 



q:P,q',k' 



X {2Trf S^{k + q~ p){2^Y5^{k' + q -p) 



1 



x//0(fc,g)/^'H-(l-e-r*) 



(11) 
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FIG. 1: 1-loop contribution to the self-energies Ep of the 
Majorana neutrino A^. 

Note that the integrand is now proportional to the aver- 
aged matrix element \M{l(l) 14>)\'^ = 2k ■ k'{X^X)ii/M'^ 
(cf. [4]), which involves the product of the 4- vectors k and 
k'. At low temperatures, T ^ M, the integrand falls off 
like e~^'^ < e~^^^ , i.e., the generated lepton asymmetry 
is strongly suppressed. 

Let us now consider spectral function and statistical 
propagator (cf. [l^) for the heavy Majorana neutrino, 

G-^ixi,X2) =i{{N^{x,),Np{x2)}) , (12) 

Gtpixi,X2) = ^{[N^{x^),Npix2)]) , (13) 

which satisfy the Kadanoff-Baym equations [1, [l3| 

C(z7°at,-p7-Af)Gp(ti-i2) = 

- / 'dt'I]p(ii -t')Gp(t'-t2) , (14) 

+ r dt'Y.+{t^-t')G-{t' 

Jo 

dt'j::p{ti~t')G+{t',t2) . (15) 



Here we have assumed spatial homogeneity and per- 
formed a Fourier transform. The 1-loop contribution to 
the self-energies Ep is shown in Fig. 1. For small cou- 
plings, A ^ 1, leading to a small width F ^ M, explicit 
solutions of the Kadanoff-Baym equations can be found 
in the Breit-Wigner approximation [lo| . 



Gp(y) 
Gt{t,y) 



i^Q cos{ujy) + 
X e-rii/l/2c-i , 
170 sm{ujy) - 
tanh (^) 



M 



P7 



sm{ojy) 



M - p7 



cos{ujy) 



(16) 



(17) 



-rt 



G- 



where t ^ {h + t2)/2, y ^ ti - t2 and F = F,3(w) 
For large t, Gp(t,y) approaches the equilibrium solution 
Gp'^'^iy), and for small temperatures, i.e., large (3, it be- 
comes the vacuum solution Gp'"°''^{y). Note that the so- 
lution p7|l for Gp{t^y) satisfies the initial condition 



G+(0,0) = G+™^(0) , 



(18) 
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FIG. 2: 2-loop contributions to the lepton self-energies Ilj^, 
which lead to non-zero lepton number densities. 



which is the analogue of the initial condition /Ar(0, ui) — 
for the distribution function. For spectral function and 
statistical propagator of lepton and Higgs fields we shall 
use the free equilibrium expressions S^{y) and Aq(y), 
respectively. The thermal equilibrium is assumed to be 
established by Standard Model interactions. 

We are now ready to calculate the lepton asymme- 
try which is generated during the approach of the right- 
handed neutrino N to equilibrium. The 'lepton number 
matrix' is obtained from the statistical propagator of the 
lepton fields, 



(19) 



One easily verifies that for free fields in equilibrium, 
Lkii\t=t' = fii{k) — fii{k), which vanishes for zero chem- 
ical potential. To leading order in A, a flavour non- 
diagonal asymmetry is generated by the 2-loop self- 
energies shown in Fig. 2 (cf. Q)- Solving the Kadanofi- 
Baym equation for to first order in the self-energy 
n^, one finds after some algebra 



-S^{t2-h)U+^^{h,h)] . (20) 

A non-zero asymmetry is generated by the departure 
of Gp from equilibrium. Due to the chiral couplings only 
a chiral projection of Gp contributes, 



Gp{t,y)PL - Pl (G+(i,y) - G+^^y)) Pl 
G^{t,y)^-cos{ujy)Q{u)e-^' . 



After a lengthy calculation one obtains for the lepton 
number matrix to leading order in the width F (To| : 



Lkij{t, t) = -tij Stt 



k-k' 
I kk'io 



{{u-k-qY + ^){{u:-k' -q'Y + ^-^) 
X fi4k,q)fi4k',q')r^^{u) 
X (cos[(fc + q-k' - q')t] + e"" (22) 

— (cos[(cj — fc - q)t] + cos[(cj - k' - q')t\)e~'^^ , 

where p = q + k = q' + k'. This expression for the 
lepton asymmetries, generated by quantum interference 
in a thermal bath, is the main result of this Letter. 

The expression (|22p contains a logarithmic divergence 
0{X'^). It has to be combined with other divergent terms 
which have been neglected in ([^^ since they do not 
contribute to leading order in F. The total divergence 
will be subtracted by a counter term. The dominant 
finite contribution to the integral stems from momenta 
k + q^k' + q'^^uj. In the zero-width limit F — 0, with 
Tt fixed, the integrand is 0{1/T). Since F oc A^, this 
dominant contribution to the integral is 0{X^). 

It is interesting to compare the diagonal elements of 
the lepton number matrix L^ijit, t) with the distribution 
functions fLiit,k) given in pTjl . As expected, the same 
CP asymmetries en appear, whereas the dependence of 
the integrands on time and temperature is different. The 
reason for the different temperature dependence is the 
fact that the matrix elements in the Boltzmann equations 
were calculated at zero temperature. Hence, the factor 
fiiji{k',q') is missing in (jlip . Since Boltzmann equations 
are local in time whereas Kadanoff-Baym equations con- 
tain 'memory effects', the different time dependence of 
the asymmetries is also expected. One consequence is 
that dtfLi{t,k)\t=o ^ 0, whereas dtL]^ij{t,t)\t=Q = 0. 
Particularly important are off-shell effects in (|22l) , which 
lead to terms oscillating in time. 

The additional thermal correction fi^{k', q') is linear in 
the distribution functions // and f^, in contrast to results 
obtained in the first two papers in f(|, but in agreement 
with 01 and [H. The memory and off-shell effects found 
in fi\ are qualitatively different from our result (|^^ . 

It is instructive to consider the approximations, even 
if they might not be justified, which lead from to 
the result of Boltzmann equations. Neglecting off-shell 
effects, i.e., imposing uj = k + q = k' + q', the cosines are 
replaced by one; performing then the zero- width approx- 
imation F ^ 0, with Tt fixed, the integral becomes 



IGtt 



k ■ k' 



/q,q',p,k' 

X {2TT)'^S\k + q- p){2TT)^S\k' + q -p) 
X .fi<p{k,q)fi<t,{k\q )fi^{uj) 



(21) 



1 _rtx2 

X - (1 -e 2 



(23) 
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Except for the factor fi^{k',q'), the only difference com- 
pared to the solution ([TT|) of the Boltzmann equations 
is the time dependence. It is obvious from Fig. 2 and 
Eq. ([^0]) that in the quantum theory the generation of 
the lepton asymmetry is nonlocal in time. This leads 
to the square of the exponential fall-off in (|23|) . On the 
contrary, in the Boltzmann equations the asymmetry is 
generated locally in time yielding a simple exponential 
behaviour. The difference can be numerically important 
at cosmologically relevant times ~ l/F. 

The calculations leading to Eq. also demonstrate 
that the result for the lepton asymmetry will be signifi- 
cantly modified by the thermal damping rates for lepton 
and Hig gs fields in the plasma. These thermal widths 
(cf. [15]) are known to be much larger than the decay 
width of the Majorana neutrino: Ti ^ g'^T > A^M 

for M < T. For quantum interferences the thermal 
damping rates are qualitatively more important than 
thermal masses which, for simplicity, we ignore in the fol- 
lowing. Including naively thermal widths for lepton and 
Higgs fields in the self-energy Ilj^ , one obtains instead of 
to leading order in these widths (Fj^ = F/ -f F^) 

jri^F^ 

{{c-k- qy + iF2^)((c. -k'- q'Y + \Tl) 
X Mk,q)fi^{k',q')f%\uj) 

X i (1 - e-^') . (24) 
Note that the factors oscillating in time have disap- 



peared. The thermal widths of lepton and Higgs fields 
have led to a behaviour which is local in time. In the 
zero-width limit one now obtains the result ()II[) of the 
Boltzmann equations except for the thermal correction 
factor fi^{k',q'). We emphasize that ([^^ is speculative 
at present, and it remains to be seen whether it follows 
from a solid calculation which includes gauge interactions 
in a systematic way. 

We have studied the generation of a lepton asymmetry 
at constant temperature as the heavy Majorana neutrino 
approaches thermal equilibrium. For cosmological lepto- 
genesis one has to calculate the lepton asymmetry in the 
case of decreasing temperature, which is caused by the 
expansion of the universe. Analyses based on Boltzmann 
equations suggest that both asymmetries are of compa- 
rable size for T ~ M 01 . A more detailed discussion will 
be given in [lo| . 

In this Letter we have compared lepton asymmetries 
calculated on the basis of Boltzmann equations with 
those obtained from Kadanoff-Baym equations. Our 
discussion illustrates, that even ignoring spectator and 
flavour effects jT^, current leptogenesis calculations 
have an uncertainty of at least one order of magnitude. 
Particularly urgent is the inclusion of gauge interac- 
tions with the thermal bath in a full quantum calculation. 
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